FIAN/TD-19/00 
ITEP-TH-66/00 
hep-th/0011222 



Seiberg-Witten Toda Chains and N=l SQCD ft 



A.Marshakov J 

• Theory Department, Lebedev Physics Institute, Moscow 117924, Russia 

f*^) | and 

ITEP, Moscow 117259, Russia, 

(N 

> 
O 

, We consider the Seiberg-Witten Toda chains arising in the context of exact solutions to Af = 2 

SUSY Yang-Mills and their relation to the properties of TV = 1 SUSY gauge theories. In particular, 
we discuss their " perturbative" and "solitonic" degenerations and demonstrate some relations of the 

■ latter ones to the confining properties of Af = 1 vacua. 

> 

(N ! 1 Introduction 

<N 

The Af = 2 SUSY gauge theories possess continuous parametric families of vacua, labeled by VEV's of complex 
scalar fields An exciting result is that the masses of BPS states and low-energy effective actions for massless 
fields in these theories can be found exactly [Q and are described in terms of integrable systems so that the 
moduli of Af — 2 theories are identified with moduli of complex structure of the corresponding spectral curves 
[2| or integrals of motion of completely integrable models || . The corresponding family of integrable models 
O in the most simple case consists of the Toda chains || of exponentially interacting particles. 

These Af = 2 SUSY gauge theories are even relatively realistic from the point of view of their physical spectra, 
containing propagating massless abelian gauge fields, but, as a consequence of extended supersymmetry they 
possess also massless scalars - the excitations of moduli, whose presence contradicts to what one should expect 
in reality. The situation is seriously improved when Af — 2 SUSY is violated down to Af = 1, then in contrast 
to Af = 2 situation, in Af = 1 case a superpotential W / ( < &) may be generated which fixes the VEV of scalars to 
the extrema dW — and one ends up with some restricted set (usually finite number) of vacua where there is 
no more place for light moduli scalar excitations. It is widely believed that Af = 1 gauge theories are already 
almost the theories close to real QCD, therefore investigation of the properties of these vacua is an important 
problem. 

Below, we will try to present some arguments that the Seiberg-Witten integrable systems, arose originally 
and associated usually to the Af = 2 situation, may also play some role in investigating properties of Af = 1 
gauge theories. We will start with generic properties of SUSY algebra in supersymmetric Yang-Mills theories 
and demonstrate that already at this level one may expect appearance of certain hamiltonian and/or integrable 
structures. Then we will pass to direct investigation of particular limits of the Seiberg-Witten integrable systems, 
since it is clear from the above reasoning that more realistic Af = 1 picture appears only after some conditions 
on moduli are imposed. We will discuss, in particular, degenerations of the smooth Seiberg-Witten curves and 
demonstrate that, in fact, these degenerations can be separated into two different classes, the properties of 
integrable systems in these limits were studied in ||. It turns out, that the structure of an integrable system 
implies a presence of a kind of " regularizing" multiplet of fundamental matter - to be identified with the set 
of Baker-Akhiezer functions of an integrable model, which can be explicitly constructed and their properties 
will be discussed below. In particular, we will see that the perturbative Baker-Akhiezer functions satisfy the 
condition, which may be thought of as a (holomorphic square root of) the D-term condition for the fundamental 
matter. We complete our discussion by few remarks about the "internal" properties of the Toda chain solutions. 

1 Based on the talks given at the E.S.Fradkin Memorial Conference, Moscow, June 2000 and NATO workshop Integrable 
Hierarchies and Modern Physics Theories, Chicago, July 2000 . 

2 e-mail address: mars@lpi.ru, mars@gate.itep.ru, marshakov@nbivms.nbi.dk 
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Figure 1: Smooth Seiberg-Witten curve for the Af = 2 SU(N) pure gauge theory associated with periodic Toda 
chain. 



2 Supersymmetry and integrable structures 

The basic understanding of integrable structure in Af = 2 SUSY gauge theories appears already at the level of the 
(extended) SUSY algebra. Indeed Af = 2 SUSY algebra in 4D is generated by the supercurrents Qa = J d 3 xJ A 
(A = 1,2) and their complex conjugated. These currents differ from each other by exchanging the fermions 
{ipA~\ an d have general structure 

J° A = (ii A F) 8 OaD) 8 e AB $ B D$) 8 e AS (V> B F) (1) 

where F = DA is gauge field strength, D Q ($ t , $) = [$t, $] Q is the D-term and F Q ($) = J|£ is the F-term^. 
From this form and canonical commutation relations it is clear that the anticommutator of Qi and Q2 would 
contain the following contributions 



{Qi,Q 2 }= {$aD$), (4>aF)} 

and these terms lead to the central charges || 



(2) 

A=l:2 



}— j D$ A {F + *F) (3) 



^SUSY — -y 

9ym 

In the case of point-like charged objects the 3-dimcnsional integral in the last formula can be rewritten using 
the Stokes theorem as an integral over sphere S 2 at spatial infinity, giving rise to Zsusy ~ o,(n e + rn m ), where 
a r~j (<j>) is VEV of scalar and r = 1? H — i— is (complexified) coupling constant. However, if we originally 

consider instead of R 3 a compactified space R 2 x S 1 , the 2-dimensional integral is rather taken over S 1 x S 1 
and can be presented as 

Zsusy = [ d 2 xe ij di<$>d 3 A = [ 5$ A 6 A (4) 

is'xS 1 J2-cycle 

The last integral may be, in fact, identified with the integral of holomorphic symplectic form <5$ A 5 A where the 
appearance of extra moduli parameters A is related to the Wilson-Polyakov loops of gauge fields, complexified 
by dual photons and becomes clear if one considers compactification of 4D theory down to 3D plus possible 
one compact dimension J7J. The loops of gauge field break the invariance in colour space and, in particular, 
this leads to replacement of the diagonal matrix of the VEV 4> by, in the simplest case of Toda chains " almost 
diagonal" , Lax operator C of an integrable model. 

Moreover, in the compactified case one may violate Af = 2 4D SUSY by the Scherk-Schwarz mechanism, 
which leads to generation of superpotentials, in spirit of M. These superpotentials exactly coincide with the 
Hamiltonians of the Seiberg-Witten integrable systems P>|9]. It means, in particular, that integrable models 
are even more visible in compactified Af = 1 theory (Af = 2 in 3D terms), than in 4D Af = 2 case. 

SUSY breaking from Af =2 down to TV = 1 (in 4D terms) is usually associated with degeneration 0, [h]. |ll| 
of smooth Seiberg-Witten curve for Af = 2, sec fig. [j] In fact not any degeneration of a curve corresponds 
to Af = 1 vacuum. Roughly speaking there are two classes of degenerations, the "open" or perturbative and 
solitonic, if speaking in the language of integrable systems, the same may be thought of as asymptotically free 
and zero-charge from the point of view of corresponding quantum field theories. Only the latter class of (solitonic 

3 The superpotential terms are, in fact, absent in Af = 2 theory, but they arise, in general, after breaking SUSY down to Af = 1, 
which is important in the discussion below. 
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Figure 2: A perturbative degeneration of the Seiberg-Witten curve from fig. shrinks into its solitonic limit. 

or zero-charge) degenerations may play the role of set of vacua for M = 1 SUSY gauge theories. Qualitatively 
the difference is clearly seen in the brane picture for the Seiberg-Witten Toda chain curves |l2|, [u| [II], [t5| , the 
brane with positive tension would like to "shrink" to the zero-charge points, even if you start from different, 
perturbative degeneration (see fig. |^). 



3 Integrable Toda chains from weak coupling 

We start from perturbative degeneration, where J\f = 2 SUSY effective actions and integrable systems behind 
are completely determined by the 1-loop contributions. In the example of pure J\f = 2 gauge theory with 
SU(N) gauge group [||, the scalar field $ acquires nonzero VEV's $ = diag(ai, . . . , ajy) extremizing the 
potential Tr[<f>, $^] 2 . This (generically) breaks the SU(N) gauge group due to the interaction term ([A^, = 
[A^iai — Oj)) 2 and the masses of lU-bosons (and their superpartners) are classically given by ciy = — cij or 

-rpcrt 



dS peit =f \dlogw=<t> \d log P N (A) (5) 

for a particular "figure-of-eight" like contour Cy around the roots A = and A = a,j of the generating 
polynomial 

N 

w = P N (\)= det N (\-<Z>) = l[(\-a i ), (6) 

i=l 

This means that the contour integrals in the complex A-plane (^) for such contours around these singular points 
(the roots of the polynomial (||)) give one set of the BPS masses in the Seiberg-Witten theory. Another set 
of masses (monopoles) may be associated to dual contours starting and ending at the points ai where dS pcrt 
has singularities, taking into account that divergences are absorbed by a renormalization of the couplings. The 
integrals over the contours on the marked plane are in fact contour integrals on a degenerate Ricmann surface 
depicted at the top of fig. ^|: if one cuts it (vertically) into two pieces, connected by very narrow tubes, one gets 
exactly two copies of the desired A-plane. 

The N = 2 perturbative effective action (or prepotential) T and the set of effective couplings —iuTjk = 

dTdat- ' are Si yen ^ the 1-loop diagrams with the result ^| 

i k 2m ^ 8 A 2 2iri 8 A 2 (7 > 

(masscsM) jfc 

4 Instanton contributions to the prepotential which are proportional to powers of A 2JV are suppressed, and one keeps only the 
terms proportional to log A or the bare coupling t. 
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where the scale parameter A = Aqcd so introduced may be related to the bare coupling r. 

The connection with integrable systems comes by interpreting the differential dS pcrt with the generating 
differential of a Hamiltonian system - an open Toda chain or Toda molecule 0. Remarkably, this interpretation 
extends beyond the perturbative regime, the "regularized" curve for (^) (see hg. [l| ) 



A2W 

w+ = Pn(X) (8) 



is the spectral curve of a periodic Toda chain (an analog of the Sine-Gordon model in the same sense as the 
Toda molecule is a multiparticle analog of the Liouville model) providing the exact non-perturbative solution 
taking into account all instantonic contributions. The residue formulas (||) for the smooth curve (||) of genus 
N — 1 turn into the action period integrals 



i t i . dw 



(a,a")=* A— (9) 



/ A,B 

so that the set of effective charges (Q) is generalized to the period matrix of (|J) 

daP d 2 T 

r« = = i,j = i,...,N-i (io) 

oaj oaiOaj 

where T now stands for non-perturbative prepotential or effective action. 

4 More on Toda: periodic chain and molecule 

The dynamics of the integrable system itself may be again seen only in the compactified version Q when the 
number of moduli of the gauge theory jumps up to the whole phase space of an integrable model. In this case 
the angle variables play the role of extra moduli related directly to the hamiltonian " times" . This relation is 
especially simple in perturbative regime, when the basis of (normalized) " holomorphic" differentials 

2iridu)i = (11) 

A - a, 

and corresponding angle variables (originally defined as co-ordinates in the Jacobian of smooth curve of 
fig. [l]) are 



N r p k ri\ 00 N r p i> j\ 

k=l JF « 1 n=0 k=l JF ° 



ATlogA + V<* r , 

A fc ^A(P )=oo ^— ; 



(12) 



for some divisor {P&} with = A(P^) and 

^E A r ( 13 ) 



N 

[ 

t n — 



n 

k=l 



are the "hamiltonian" times in Miwa variables. It means, that after renormalization z% — > Zi — j^j log A one gets 

2iriz l — a t t + ^ a"t„ 

n>l 

A natural ingredient of the integrable formulation is the set of the Baker- Akhiezer functions {^i} (and {'I'}), 
which can be thought of as introducing of "regularizing" fundamental matter. Like quark multiplets in SUSY 
gauge theories the Baker- Akhiezer functions always appear "in pairs" and carry a single index of lowest fun- 
damental, N - for {^i}, or N - for {^i}, representation of gauge group. Again, in 4D theory they are, up to 

5 From the simplest SU{2) pure gauge theory [0, with the curve w = A 2 — h, h = ^Tr"!? 2 = a? it becomes clear, that the 
masses feb are defined by the contour integrals of dS pcIt = ^ d \_ -| — >"!■>■ w j 1 j c j 1 can b e interpreted as an integration of open Toda- 

l-J \ — y/h \-\-y/h 

Liouville model with Hamiltonian h = p 2 4- e 2q , provided A — > p, w — > e 2q . The same is true in general, that is, integration of the 
canonical differential dS = Ipdq over the trajectories of the Toda chain solutions gives rise, for various (complexified) trajectories, 
to the BPS masses in the Seiberg-Witten theory |3| . 
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gauge transformation, just the eigenvectors of diagonal matrix but in compactified version, when $ should 
be replaced by the Lax operator C, they become nontrivial. 

On degenerate spectral curve (||) the set of the Baker-Akhiezer functions can be defined by the following 
analytic requirements ||: 

• This is the set of functions = ^(A), k = 0, . . . , N — 1 which have exactly k zeroes on rational curve 
(JsJ) and a single pole of order k at A = oo. For k > N they can be defined by ^k+N = w^k- 

• These functions can be constructed as linear combinations of Yi% 1 < <i k (^ ~ a ii) = Ilje/C^ — a i) or > 
vice versa, certain linear combinations of vanish exactly at the chosen sets I oi k points a^, . . . ,cti h , 
*!<...< ife. One can see that the Baker-Akhiezer functions, indeed, satisfy for e Zi = e 2 '* %Zi X\j^ i -^— 

N 

E^K) = (-) k+1 J2 Re*' II a ^ = ( 15 ) 

j=i \i\=k+nei ijei 

It is easy to see, that so defined Baker-Akhiezer functions are related to the tau-functions and can be written 
(up to normalization) as 

**(A) = a* e k L - E Jr) = n^CA - ai )H±- (is) 

V n J iei jel 13 

where we write all a%j in such order that i < j. The sum ^ fe is taken over all sets of |7| = k indices with I 
denoting the complementary set to /. The tau-functions 

~. X 2tti V W 2, — Z7r (ni—nj) 2 T n . n . 

=E fc iP^n^ 

iei jel 

are nothing but degenerate theta-functions for the Toda molecule system where in the perturbative limit A — > 
the period matrix is given by (Q) with the terms proportional to log A absorbed into some renormalizations: 
Zi — ► Zi — log A, Ofc — * A fe Ofe Q They can be obtained as degenerations of the theta-functions on (Jacobian 
of) smooth curve (^) and their properties can be found, for example, in The formulas ( [l6| ) can be written 
down even more explicitly, say for the N = 2 case one gets 

•.m-£(-(-e£)+-4-?h-£))- 

p z -X- p z 

(A — atanhz) = e~ 9 (A — p) 



a 

providing the relations Oi = \{e z + e~ z ) = e~ q and p = atanhz ^. One may also easily check that the vector 

*(A) = ( 1 ) (19) 
V *i(A) J 

is an eigenvector of the open Toda Lax operator C^f = A\t 

( p e« \ 

C_ U -J (20) 

provided by A 2 = p 2 + e 2g = h = a 2 , which is related to the VEV of 4D scalar <£> by C — U^U^ 1 , where matrix 
U is constructed from the Baker- Akhiczcr functions. 



6 The TV-dimensional vectors like {'I'fc} are labeled here for simplicity of notations as fe = 0, . . . , N — 1, unlike usual conventions 
in group theory. 

7 The last renormalization could be understood as a logarithmic renormalization of the bare coupling r. 

8 Note, that "3/2 = A 2 — a 2 = when A = ±a belongs for a spectrum of the Toda molecule. This is a common feature for all 



*jv=nHi( A -^)- 
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The dual Baker- Akhiczcr functions {'I'fc(A)} (which may be thought as coming from the "dual" sheet of 
degenerate Toda chain curve Jl6[) can be written, up to normalization, again in terms of the theta- functions as 

\ n / ieJ jg/ J 

For example 

1 



riili( A - a < 



W — ■ 1 _ Ef^e 2 — U'j^ (22) 



n l =i( A - a * 



so that for the N = 2 case one gets 



A 2 - a 2 



*o(A) = A" 1 - (e 2+ £-> * + e- 2+ S„H"^) = A ' la + = (23) 
a \ J a 

'e z +e~ z \ A + atanhz A+p 



A 2 — a 2 A 2 — a 2 

It is easy to see, for ([L8|) and (|23| ) in particular, that one may indeed impose gluing conditions 

V i (a)=CV i {a) (24) 

for some coefficient C. In our basis ([H]) for the dual Baker- Akhiezer functions this coefficient has to be 
nontrivial, since ^-(a) are naively singular exactly in the points of the spectrum, but adjusting C, or changing 
the normalization of (|2l|) we see, that the gluing condition ( p4] ) indeed holds. 

These are the gluing conditions at weak coupling - and they may be thought of as taking holomorphic square 
root from the D-term condition for a single flavour of "regularizing" fundamental matter. 



5 The strong coupling limit 

The quantum moduli space of the 4D pure SU(N) M = 2 SYM has N maximally singular points at which N— 1 
monopoles become simultaneously massless - the confining vacua of an M = 1 theory g, |, ||| , where the 
dual variables af = are appropriate to describe the prepotential. The J\f = 2 spectral curve (||) at these 
points is described in terms of a Chebyshev polynomial (in this section we put A = 1) 

p Choby S hcv (A) = 2cQshz (25) 

with a solution A = 2 cosh (j^ + i^ff-) (k = 0, . . . N — 1). These are the N = 1 points of the theory, related by 
a Zjv symmetry. The hyperelliptic form of the spectral curve (j^) (recall y — w — l/w) is now 

N-l 

y 2 = P N (X) 2 - 4 = (A 2 - 4) [] (A - 2 cos ^) 2 (26) 

j=i 

which is a "solitonic" curve in the periodic Toda chain now the corresponding B periods have collapsed, 

in contrast to the perturbative case. Introducing 

A = 2cosh^-=£ + r 1 - (27) 



one maps the 2-sheeted cover of the A-plane y = \/A 2 — 4 to a cylinder with co-ordinate £. Thus eqs. (|25|), ( p6| ) 
describe analytically a cylinder with iV — 1 distinguished pairs of points - an inside-out S-duality transform of 
the bottom curve from fig. @. It is worth noting that in contrast to perturbative regime (the upper curve on 



fig. |^) now a curve does not separate into two pieces. The canonical basis of " holomorphic" differentials in this 
limit can be can be chosen as (j = 1, . . . , N — 1) 



, D sm^f i c-g 

du)~ = —-, r~/ r = ^ — dlog —. (oe\ 

J IT r r —*Si\ , 27T7 & F-F7 



These differentials are normalized to the B-cycles, here the cycles around the marked points 

D _ X j, ,D 1 j, ,D 



dujj = f duj = - f dulj = Sij, (29) 

while certain of the A-periods (§ A dui® — f i dujf) diverge, the others (j ^ k) being given by [] 

3 ?j 



,, , ,, 1 sin 2 — k) 1 sin — k I 

TR= i dcj£ = — log 2 2 ^ -f = — log . ™\. -i. 30 

' " 2m & sin 2 ^(j + k) in B sm^(j + fc) 

This is an exact formula for the set of effective couplings in strong-coupling regime, which has "typical solitonic" 
form. The appearance of solitons at strong coupling is a sort of "intuitive explanation" of the absence of the 
point with ($) = since the solitonic points correspond, on one hand, to maximal degeneracy and, on the other 
hand, to minimum of superpotential, being an analog of the "BPS" first-order solutions in Seiberg-Witten Toda 
chains. 

The angles at strong coupling are given by 



N—i -, N-l 

n 

z 



-, iv — 1 £ ft -, oo 

where t n = ^2 k depending on whether k = 1, . . . , N — 1 go to or to oo. We observe also that the 



relation ([31]) coincides with the vacuum value of the string tension, or monopole condensate, proportional to the 
SUSY breaking parameters tfc |l0|, |l5|, [l8|, [l9| . The detailed study of the properties of these strings are beyond 
the scope of this note, a discussion of related questions can be found in po| . 

In contrast to the weak-coupling regime, the gluing conditions at the strong coupling are |L8[ |) 

*»(#) = *„($-) j = l,...,N-l. (32) 

which means that the Baker- Akhiezcr function on a cylinder remembers that it originally came from a genus 
N — 1 Riemann surface, and each pair of points = corresponds to a degenerate handle. Note, that 

the form of condition (^||) is different from ( pi]) (and in fact, is even more similar to (|l5|)) and relates the 
same Baker- Akhiezer functions in different points of the strong-coupling spectrum. There is no natural relation 
between VP and in solitonic phase, which should be probably related to the properties of fundamental matter 
multiplets at strong coupling. We are going to return to this problem elsewhere. 



6 More remarks on the Toda tau functions 

In this section we will try to make few more remarks concerning the degenerate limits of the Toda chain, which 
do not have at the moment direct relation to the properties of Af = 1 SQCD, nevertheless being of some interest 
by themselves. We believe, however, that their physics origin is related to the bulk/boundary duality between 
SUSY gauge theories and gravity. 

Commutativity of the theta-functions. 

The theta-functions in ( p"7| ) may be obtained as degenerate limit of the coefficients Ofc in the expansion 
of the genus N theta function 

e N = J2 e 2 ^ z e fe (z|T) 

9 This formula was independently proposed in JlT^. 
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We may already deduce an interesting consequence: the ratios of the coefficients 0fe = 0fc(z|T) Poisson commute 
(see 1 22 and references therein), 

|,|£} =0 , V,,,,,',/. ,34) 
The Poisson bracket here is that corresponding to the symplectic form 

N-l N-l 

n Toda = ^ dz, A dat = dqi A d Pl , (35) 

i=l i=l 

where qi and pi are co-ordinates and momenta of the Toda chain and we are working in the centre of mass frame. 
The Poisson commutativity of the ratios (|34|) follows from the solution of the periodic Toda chain p3l p3 



6; 



k=j+l 



,1k 



(36) 



(Here Qq — &n and the 0fc are the Toda chain tau-functions depending on the discrete time k, the number 
of the particle.) Now because the coordinates qi obviously Poisson commute, {q%,qj} — 0, we deduce (|34|). 
This expression gives a precise formulation of the old expectation that the Toda chain tau-functions Poisson 
commute with each other. In the case of periodic Toda chain, however, only the ratios are commutative (theta- 
functions themselves, co-ordinates in projective space, arc defined only up to common multiplication), but in 
the degenerate perturbative limit Go — 0^ = 1 and, thus, tau-functions commute themselves. Moreover, in 
the perturbative limit formulas (17) and (|36j) give an explicit solution to the Toda molecule Q in terms of the 
action-angle variables. 

Relation to matrix models. 

It is widely known now that the partition functions of the matrix models are tau functions of Toda chains 
I PS) . Expression (|l7|), being a finite sum, still is a sort of partition function, when the sum is taken over the 
discrete spectrum of the finite Toda molecule. It is easy to demonstrate, that, indeed this tau function can be 
considered as a " discretized" analog of the partition function of hcrmitian 1-matrix model. 



For example, one can verify that 



where the "moment matrix" K„ 



Kn 



e fe = det K n+m 

X n,m— 

K n+m is "an average" w.r.t. 0i, i.e. 

N 
i=l 



,Zi 



(37) 



(38) 



(cf. with ^||, |2(|). The simplest proof Q is based on the Cauchy-Binet formula for the k x N (k < N) 
rectangular matrices A n i and Sj m (i = 1, . . . , N, n,m — 1, . . . , k) 



N 



n=l,. 



Edet A 
fcxfe 

n,m=l,...,k h<-..<ik 

so that 

Qk = 

^ i -i. 

iei jei J iei i^jei 

is indeed a "discretized" analogue of the tau-function of 1-matrix model J 



det Bj 



kxk 



det V A ni B, 

kxk \ ^ 

\i=l 



iei; n=l, 



e 



n 



i j 



J2\I\=k TliEl Ci llj<j; i,j£l a ij — h<...<i k C*l 

= J2t 1 <...<i k Cii ■■ ■ C*k detfexfe a" | ig/ . n=1 ^ k 
= detfe X fe (YZi C*K +m - 2 ' 



i£l i<j; i,j€l 

since for any coefficients Ci 

\'2 



det 



r, 



(a in ~ di m ) 



,m — 1 1 



(39) 



(40) 



(41) 



and, substituting Ci = e z \ one arrives at eqs. 



n,m— 1,. . . ,k 



E 



o k = > 

/: i i < . . . < i fe 



e Z ia n+m-2 



det K n+m . 

kxk 



(42) 



The " non-discretized" (though still " discrete" ) 1-matrix model should be restored in the N — > oo limit of the 
Toda molecule, under condition that the spectrum {a.;} fills in the whole real axis. Then instead of the Cauchy- 
Binet formula (rectangular matrices become infinite in one of the direction) one can effectively use the technique 
of the orthogonal polynomials p5], |26]| . 
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7 Conclusion 



We have discussed just a brief list of questions related to the problems of AT = 1 vacua is SUSY gauge theories. 
However, already at this level we can see that there are remarkable coincidence between the properties of J\f = 1 
theories and Seiberg-Witten Toda chains " at strong coupling" . They include the relations between the Baker- 
Akhiezer functions and the formulas for the phases of solitons, which look like being "in charge" for confining 
strings in SUSY gauge theories. It is necessary also to point out that perturbative integrable system is "stable" 
under small changes of moduli a i; while in strong coupling regime all moduli are completely "confined" to 
certain values. 
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